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ABSTRACT 

There is an apparent dichotomy between the metal-poor ([Fe/H] ^ —2) yet carbon- 
normal giants and their carbon- rich counterparts. The former undergo significant de- 
pletion of carbon on the red giant branch after they have undergone first dredge-up, 
whereas the latter do not appear to experience significant depletion. We investigate 
this in the context that the extra mixing occurs via the thermohaline instability that 
arises due to the burning of '^He. We present the evolution of [C/Fe], [N/Fe] and 
12q^13q ^]-^j.gg models: a carbon-normal metal-poor star, and two stars that have 
accreted material from a 1.5 M© AGB companion, one having received 0.01 M© of 
material and the other having received 0.1 M©. We find the behaviour of the carbon- 
normal metal-poor stars is well reproduced by this mechanism. In addition, our models 
also show that the efficiency of carbon-depletion is significantly reduced in carbon-rich 
stars. This extra-mixing mechanism is able to reproduce the observed properties of 
both carbon-normal and carbon-rich stars. 
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1 INTRODUCTION 

The pursuit of the most metal-poor stars in the Universe 
is one of the most active fields of research in modern stel- 
lar astronomy. Several large-scale surveys, such as the 'First 
Stars ' program of Cayrel and collaborators (e.g. lCavrel et all 
|2004| ). the Hamburg/ESO survey of Christlieb and collabo- 
rators fe.g. lChristliebll2003t) an d the HK survey of Beers and 
collaborators (e.g. Beers et al.lll99j ) have yielded a wealth 
of abundance determinations. This accumulation of data 
provides us with statistically meaningful populations that 
can be compared to theoretical predictions. One result of 
such survey work is that two populations of stars can be 
defined based on their carbon abundances. There is a dis- 
tinct separation between those stars that are carbon-rich 
and those that are carbon- normal. The divi ding line is drawn 
at [C/Fej3>-f 1.0 l|Beers fc Christfiebll2005D and the carbon- 
rich population accou nts for a substantial f raction of all the 
metal-poor stars (e.g. iLucatello et al.ll200^ ). 

There is an interesting dichotomy between the be- 
haviour of carbon-normal and carbon-enhanced metal-poor 
(CEMP) stars. This is shown in Figure [1] which shows mea- 
surements of [C/Fe] as a function of surface gravity, logjg g. 
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Stars have higher surface gravities when they are more com- 
pact. On the main sequence logj^g 5 ~ 5 and as the star 
expands, logj^y g drops. By the time logj^Q g « 3, the star 
has reached the giant branch and its convective envelope 
deepens, bringing CN-cycled material to the surface. This is 
first dredge-up. The star continues to expand and by the tip 
of the red giant branch, logjg g « 0.5. The carbon-normal 
stars, once they reach a surface gravity of around logjg g ~ 2 
(i.e. stars that have passed through first dredge-up and are 
on the upper part of the red giant branch), show a marked 
depletion in their surface carbon abundance. [C/Fe] drops 
by nearly one dex between logjQ g » 2 — 0.5. On the ba sis of 
this drop in the [C/Fe] abundance. lAoki et al.l (|2007l ) sug- 
gested that the definition of carbon-enhancement should be 
revised to take account of this depletion. They suggest that 
the appropriate criteria should be: 

(i) [C/Fe] ^ -fO.7, for stars with log (L/Lq) sC 2.3 
(u) [C/Fe] ^ +3.0 - log (L/Lo), for stars with 
log (L/Lo) > 2.3 

While there is clear evidence for the depletion of car- 
bon in the C-normal population, the carbon-rich metal-poor 
stars show no sig n of carbon depletion when logm.'? < 2. 
This was noted bv lDenissenkov fc Pinsonneaultl l|2008bl ). In 
fact, there is a noticeable dearth of stars in the triangu- 
lar regions whose apexes are at (logj^Q g, [C/Fe]) = (2,1), 
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Figure 1. [C/Fe] as a function of log^n fl for meta l-poor stars 
extracted from the SAGA database IjSuda et al.ll2008h . Light grey 
crosses represent stars that only have [C/Fe] measured, while dark 
open triangles represent stars that have measurements of [C/Fe] 
and [N/Fe], The black open squares are those stars for which 
we have data for [C/Fe], [N/Fe] and i^c/Wc. The solid line at 
[C/Fe] = 1 denotes the dividing line between carbon- normal and 
carbon-rich metal-poor stars. 



(0.5,1), (0.5,-0.5). This is unlikely to be a selection effect as 
observations are sensitive to this region of parameter space 
iLai et all [20071 ). It therefore seems that carbon-enhanced 
metal-poor stars do not suffer the same degree of carbon- 
depletion on the upper part of the red giant branch as their 
carbon-normal cousins do. 

The existence of extra mixing i n red giant bran c h star s 
is observationally w ell established. IPearborn et al.l (|l975l) , 
iTomkin etaP (|l976l l and more recently [Keller et all (|200lh 
used "C/"C ratios as an indicator of the efficiency of mix- 
ing associated with first dredge up. The discrepancy between 
canonical models and observations were the first indications 
that sta rs might undergo some fo rm o f extra mixing. Furthe r 
work bv lCharbonnel et all (|l998l ) and lSmilianic et al.l (|2008l ) 
showed that there are unexplained increases in N whilst ^Li, 
C and ^^C/^'^C all decrease. Systematic studies of various 
stellar environments - including field stars, globular clus- 
ters, open clusters and the Large Magellanic Cloud - have 
shown evidence of extra mixing operat ing universally across 
a range of masses and metallicities (see lSmilianic et al ]|2008l . 
and references therein). T he efficiency of th e mixing has 
also been well studied with lKeller et all (l200ll) determining 
mixin g is more efficient at lower metallicities. iMartell et al.l 
(|2008t ) show progressive C depletion with increasing lumi- 
nosity along the giant branch suggesting that the mixing 
occurs throughout a star's ascent of the giant branch. 

The mechanism for this additional mixing has proved 
somewha t elusive. E xtra mixing caused by rotation was pro- 

f osed by IPaczyiiskil (Il97 3'l and again by CottrcU & N orrisI 
19781 ') to explain low C/N ratios in young sta r s and weak 
G band stars respectively. ISweieart fc Mengell ()l979l ) were 
the first to investigate the role of meridional circulation on 
the RGB and this has been an avenue o f pursuit for many 
authors si nce. iDenissenkov fc ToiI3 (|2000l ) , applying the for- 
malism of iMaeder fc Zahn ( 19981 ) in a post-processing ap- 
proach, showed that rotational mixing could produce abun- 
dance variations on the RGB. However, their models also 
produced anticorrelations in 0-Na and Mg-Al. It is now 



thought that these features are not caused by extra mixing 
on the giant branch, as they have been observed in turn- 
off stars in globular clusters. IDenissenkov fc VandenBerd 
(2003) subsequently improved on this work, showing that 
the abundance anomalies could be reproduced if the diffu- 
sion coefficient from their rotating model wa s multiplied by 
a fact or of 7. However, the recent models bv lPalacios et all 
(200^) cast serious doubt on the assumption that rotation 
is the mechanism responsible for extra mixing on the RGB. 
Using models in which the transport of angular momentum 
by meridional circulation and shear turbulence, in addition 
to the associated chemical mixing these processes cause, was 
self-consistently treated (i.e. the feedback between rotational 
processes and stellar structure was taken into account), these 
authors showed that rotational mixing was unable to ac- 
count for the o bserved abundance changes (see figure 15 in 
Palacios et~ail [20061 . 

Recently, an alternative mechanism has received much 
attention as the possible cause of mixing on the RGB. 
The direct numerical simulations of lEggleton et all (|2006| ) 
showed that the reaction ^He(^He,2p)*He (this reaction and 
its effect upon thermohaline mi xing in pre-m ain sequence gi- 
ants was previously studied bv lUlrich[|l972l ) could lead to a 
lowering of the mean molecular weight above the H-burning 
shell and that this would lead to mixing. In follow-up work, 
these authors modelled this process via a diffusion coeffi- 
cient which gave diffusive velocities in the range expected 
from simple physical arguments. They then went on to show 
how this could affect the isotopic ratios of various elements 
(jEggleton et al ]|2008l ). As a star evolves up the giant branch, 
its convective envelope deepens and material that has under- 
gone CN-cycling is brought to the surface. This is the pro- 
cess known as first dredge-up. First dredge-up homogenises 
the stellar envelope, resulting in a uniform mean molecu- 
lar weight /i. As the hydrogen burning shell reaches this 
homogenous region '^He starts to burn, lowering the mean 
molecular weight in this region. The resultant /i-profile is 
unstable to thermohaline mixin£[ and the material begins 
to be circulated between the burning region and the con- 
vective envelope. This connection of the burning region to 
the envelope (and hence the stellar surface) causes the sur- 
face abundances to change, with carbon becoming depleted 
and nitrogen enhanced. One particularly appealing aspect 
of this mechanism is that it is a direct consequence of the 
physics of nuclear burning and does not require a particu- 
lar rotation rate or rotational profile (or even magnetic field 
strength, if such physics were to be included) - hence the 
reason that Eg gleton and collaborat ors dubbed the process 
"compulsory" (jEggleton et al.|[20o3 ). 

The effect of this mixin g' has been investi- 

gated for a range of rn etallicities ([Eggleton et al.l [20081 : 
[Charbonnel fc Zahn|[2007l ). with the general conclusion that 
mixing is more efficient at lower metallicity, in line with ob- 
servations. In this paper, we apply this mechanism to both 
the carbon-rich and the carbon-normal metal-poor popula- 
tions to see what effects it may have. 



^ We use the terms thermohaline mixing and 5)1 mixing inter- 
changeably throughout this paper. 
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2 THE STELLAR EVOLUTION CODE 

Calculations in this work have been carried out using a mod- 
ified versio n of the STARS ste llar evolution code originally de- 
velo ped by Eggletonl l|l97lh and updated by many authors 

a). 



(e.g. lPols etalll 19951 ). The code solves the equations of stel- 



lar structure and chemical evolution in a fully simultaneous 
manner, iterating on all variables at the same time in order 
to converge a model fsee IStancliffdlioO^ . for a detailed dis- 
cussion). The version used he r e inc ludes the nucleosynthe- 
sis routines of IStancliffe et al.l (|2005) and Sta nchffe {20oE[) . 
which follow the nucleosynthesis of 40 isotopes from D to 
^^S and important iron group elements. Th e code uses the 
opacity routines of lEldridge fc ToutI (12004), which employ 
interpolation in the OPAL tables (* Iglesias fc Rogers! 119961 ) 
and which account for the variation in opacity as the C and 
O content of the material varies. In ad dition, molecular opac- 
ities a re accounted for as described bv lStancliffe fc Glebbeekl 
l|2008t ). We note that metal-poor stars are typically a- 
enhanced and so the use of opacity tables computed with 
such a composition would be appropriate. However, tables 
with both Q-enhancement and variable C and O content do 
not currently exist. 

We have evolved three models, all with a metallicity 
of Z = 10~* ([Fe/H] « —2.3): one is a standard carbon- 
normal, metal-poor model with solar-scaled abundances 
(jAnders fc Grevessd Il989l ) w hile the other two ar e post- 
accretion models taken from IStancliffe fc Glebbeekl (|2008l ) 
and Stancliffc (2009). The former has been evolved from the 
pre-main sequence and has a mass of 0.8 M©. The latter two 
have post-accretion masses of 0.8 M©, after having accreted 
0.1 and 0.01 M© of material (at a rate of lO~®M0yr"^ which 
is roughly equivalent to the rate of accretion expected in 
a wind mass transfer scenario) from a 1.5 Mq AGB com- 
panion. These two models were evolved using thermoha- 
line mixing and as such the accreted material has been 
mixed to equilibrium, with accreted material reaching a 
depth of 0.6 a nd 0.33 Mq from the surface respectively (see 
IStancliffe fc Glcbbc ck 2008 , for further details of these mod- 
els). The mixing reaches equilibrium before the end of the 
main sequence. In each case, the models have been remeshed 
to 999 meshpoints immediately after first dredge-up. Ther- 
mohaline mixing has been included via the prescription of 
iKippenhahn. Ruschenplatt fc Thomas! (|198G! ). with the dif- 
fusion coefficient multiplie d by a factor of 1000, following 
ICharbonnel fc Zah n' ('2007^. Note that this is not consistent 
with th e models of Stancliffc fc Glcbbcck (|2008l ). who use 
just the IKippenhahn et al.l (I198Q ) prescription for the mix- 
ing of accreted material. This is unlikely to effect the ex- 
tent to which accreted material is mixed during the main 
sequence as the timescale for this processes is already sig- 
nificantly faster than the nuclear timescale that governs the 
star's evolution at this point. 



3 RESULTS 

Before embarking on a description of the results of our sim- 
ulations it is necessary to discuss the sample of observations 
against which we are going to compare our models. We have 
obtained our sample from the Stell ar Abundances fo r Galac- 
tic Archaeology (SAGA) database (|Suda et al.ll2008l ). select- 
ing those stars which have [Fc/H]^ —2 and for which [C/Fe] 
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Figure 2. The evolution of the [C/Fe] abundance as a function 
of logj^Q g for three models. The solid line is for a standard metal- 
poor, carbon-normal model. The other two models have accreted 
material from a 1.5 Mq companion and this material has been 
allowed to mix via the thermohaline instability. The dotted line 
is for the case that 0.01 Mq of material has been accreted, while 
the dashed line is for the case that 0.1 Mq has been accreted. 
The initial drop in [C/Fe] at around logj^Q g x 4.5 is caused by 
thermohaline mixing acting on the accreted material. Crosses rep- 
resent objects for which we only have [C/Fe], while open triangles 
represent objects for which we have [C/Fe] and [N/Fe]. Objects 
for which [C/Fe], [N/Fe] and ^^C/^^C are known are denoted by 
open squares. 



has been measured. This gives us a sample of 621 stars. Not 
all these stars have measurements of [N/Fe] and ^^C/^'^C, so 
we subdivide our sample into 3 sets: those stars for which we 
only have a [C/Fe] determination (444 stars), those stars for 
which we have both [C/Fe] and [N/Fe] but no ^^C/^^C (122 
stars) and finally those stars for which we have a determi- 
nation for [C/Fe], [N/Fe] and ^^C/^^C (55 stars). One note 
of caution should be added: the stars we have selected come 
from a range of different surveys and so do not represent 
a homogenous sample. This is forced on us of necessity - 
we are collecting together data that was originally obtained 
for other purposes. There does not exist (to the knowledge 
of the authors) a single, homogenous set of data that mea- 
sures all the abundances we require and yet covers both the 
carbon-normal and the carbon-enhanced metal-poor stars. 

The results of the model runs are displayed in Figure (2] 
The standard carbon-normal model agrees very well with 
the data for stars with [C/Fe] < 1. Its [C/Fe] abundance 
falls by about 0.73 dex due to the effect of Sfj. mixing and 
this process sets in at log^Q g « 1.5, just as is observed in the 
data. The model which has accreted 0.01 Mq shows a much 
shallower drop in its carbon abundance with a change of 
about 0.4 dex. The model which has accreted 0.1 Mq shows 
almost no change (less than 0.08 dex) in [C/Fe] as it evolves 
to lower surface gravities. 

The variation in behaviour can be explained by look- 
ing at the evolution of the interiors of these models, par- 
ticularly at the temperature, fi and ^He abundance profiles. 
These are displayed in Figures |31 13] and [S] In each case, a 
similar reservoir of ^He is available and the dip in fi caused 
by ^He burning is roughly the same. Comparing the carbon- 
normal model (Figure (21 with the 0.01 Mq accretion case 
(Figure |4}, we see that the two models have almost identical 
temperature profiles with the /i dip occurring at a mass co- 
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Figure 3. Profiles as a function of mass for the standard carbon- 
normal model. The panels display the temperature profile (top 
panel), '^He abundance (middle panel) and fi (bottom panel) as 
the model evolves. The solid line is the earliest model, represent- 
ing the point at which 5fi mixing begins to affect the surface 
abundances. 
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Figure 4. Profiles as a function of mass for the carbon-rich model 
in which 0.01 Mq of material was accreted. The panels display the 
temperature profile (top panel), ^He abundance (middle panel) 
and /i (bottom panel) as the model evolves. The solid line is the 
earliest model, representing the point at which 5fi mixing begins 
to affect the surface abundances. 



ordinate of about 0.36 M©. However, the 0.01 M© accretion 
case evolves to the tip of the giant branch faster than the 
carbon-normal model. It behaves like a star of higher metal- 
licity and has a smaller core mass at the point of helium 
ignition. Evolution from the onset oi Sfi mixing to the tip 
of the giant branch takes 1.75 x 10^ yr in the carbon- normal 
model, while in the case of the 0.01 M0 accretion model, this 
time is just 1.4 x 10^ yr. Carbon depletion is more efficient 
at later times when the temperature at the base of the mix- 
ing region is higher, so although the time spent undergoing 
mixing differs by 20 per cent the difference in the level of 
carbon depletion is much greater, by a factor of about 2. 

One might assume that this argument would extend to 
the model in which 0.1 Mq of material has been accreted. It 
has an an even larger carbon abundance in its envelope and 
so one might expect the core to grow faster and the star to 
spend even less time undergoing Sjj. mixing. However, this 
is not the case: the evolution for the onset of Sfj, mixing 
to the tip of the giant branch takes 2.4 x lO'' yr. We must 
look more closely at the structure of this model if we are 
to understand why this model suffers less carbon depletion. 



Firstly, we note that the core mass at which mixing com- 
mences in this model is much smaller than in the previous 
two cases. In the 0.1 M© accretion model, mixing begins at 
a core mass of about 0.3 Mq, compared with 0.35 M0 in the 
other two cases. This then affects the temperature structure: 
it is significantly cooler in the 0.1 M© case (see Figure [5]). 
The reduced temperature means that less CN cycling occurs 
in the layers of the star that undergo mixing and hence there 
is less of a drop in the carbon abundance. 

Figure in] shows the evolution of [N/Fe] for our models 
alongside the [N/Fe] measurements f or our sample of stars 
extracted from the SAGA database jSuda e t al. 2008). We 
note that the 0.1 Mq model, which underwent extensive mix- 
ing on the main sequence, shows an extremely large increase 
in its nitrogen abundance at first dredge-up (logjQ (? ^ 3). 
This is because the material that was accreted has been 
mixed to a depth of around 0.6 Mq from the surface. At this 
depth the temperature is high enough for CN-cycling to oc- 
cur, co nverting the accreted C into N during the main se- 
quence (|Stancliffe et aLlboOTi ) . This material is then brought 
to the surface at first dredge-up. Each of the model se- 
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Figure 5. Profiles as a function of mass for the carbon-rich model 
in which 0.1 Mq of material was accreted. The panels display the 
temperature profile (top panel), ^He abundance (middle panel) 
and 11 (bottom panel) as the model evolves. The solid line is the 
earliest model, representing the point at which Sfi mixing begins 
to affect the surface abundances. 
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Figure 6. The evolution of the [N/Fe] abundance as a function 
of logj^Q g for three models. The solid line is for a standard metal- 
poor, carbon normal model. The other two models have accreted 
material from a 1.5 Mq companion and this material has been 
allowed to mix via the thermohaline instability. The dotted line 
is for the case that 0.01 Mq of material has been accreted, while 
the dashed line is for the case that 0.1 Mq has been accreted. Note 
that the drop in abundance at logj^Q g si 4.5 occurs as the accreted 
material is mixed into the star via thermohaline mixing. Triangles 
represent those objects for which we only have [C/Fe] and [N/Fe], 
while squares represent those objects for which we also have the 
12(~<^13q fatio. Open symbols are for C-normal objects and filled 
symbols denote C-rich objects. 
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quences then shows a rise in the nitrogen abundance when 
the Sfj, mixing occurs. The final nitrogen abundance is pro- 
portional to the quantity of carbon available and hence the 
carbon-normal model ends up being the most nitrogen poor. 
We note that the models do not cover the full spread in nitro- 
gen abundances, particularly at high surface gravities (i.e. 
on the main-sequence and at turn-ofT). However, the produc- 
tion of nitrogen from carbon does not exceed the observed 
values in any of our models. 

In Figure [7]we show the evolution of the ^^C/^^C ratio 
as a function of surface gravity for our models. The models 
that have accreted material from an AGB companion ini- 
tially have extremely high ^^C/^^C ratios (of several thou- 
sand). At first dredge-up the ratio drops dramatically in 
the 0.1 Mq accretion case because its accreted material has 
been mixed deeply into the star by thermohaline mixing 
and CN-cycled before being brought to the surface again. 
The 0.01 Mq model shows a smaller drop in its ^^C/"C 
ratio because its material is not mixed as deeply (i.e. to 
those points in the star where CN-cycling is possible) and 
the material suffers only dilution during first dredge-up. The 



Figure 7. Evolution of the ^^C/^^C ratio as a function of log^Q g 
in each of our models. The solid line indicates the carbon- normal 
model, the dotted line is the model which has accreted 0.01 Mq 
of material and the dashed line is the model which has accreted 
0.1 Mq of material. The dot-dashed line is for a model accreting 
0.1 Mq of material that has a l^c/lSC ratio of around 100. The 
open squares represent C-normal stars, while the filled squares 
are for C-rich stars. 



12q^13q ratio is substantially affected by the extra mixing 
on the upper part of the giant branch. In all cases, we find 
that this ratio is reduced to less than 10 in all our models by 
the time the top of the red giant branch is reached. This is 
unsuprising, given the quantity of nitrogen that is produced 
in these models. 

For the carbon-normal metal-poor stars, our model 
i^C/^C r atios are in reason able agreement with the obser- 
vations of ISpite et al.l l|2006l ). which are shown as crosses in 
Figure [71 The post-first dredge-up abundances are towards 
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Figure 8. Profiles of the mass fractions of the dominant CNO 
isotopes and ^He as a function of radius for our 0.8 Mq carbon- 
normal model, just prior to the onset of the extra mixing. The 
vertical line marks the radius at which the /^-inversion develops. 

the upper end of the observed ratios, but the dechne owing to 
the extra mixing is well reproduced. For the CEMP stars, we 
note o bservations of the ^^C/^'^C ratio tend t o be around 10 
or less (jSivarani et al.ll2006l : feeers et al.ll2007^ ■ This is found 
across the whole range of surface gravity, i.e. in stars from 
the main sequence to the giant branch, and hence cannot be 
held as evidence in favour of giant branch mixing in CEMP 
stars. For those CEMP stars formed by mass transfer from 
an AGB companion, there would need to be an extra mixing 
mechanism in operation during the AGB in order to prevent 
the ^^C/^'^C ratios reaching the extremely high values pre- 
dicted by our models. If a companion star accreted material 
of a much lower ^^C/^'^C ratio, it would evolve in a similar 
way to the models already presented. The dot-dashed line 
of Figure [7] shows what would happen to our model which 
accretes 0.1 M© of material if the ^^C/^'^C ratio of the ejecta 
were closer to the upper limit of the observations. The model 
suffers less of a drop in this ratio at first dredge-up because 
the envelope is closer to the equilibrium value expect from 
CN-cycling. The final ^^C/"C ratio reached at the tip of 
the giant branch is almost identical to that of the original 
model. 

We note that our results a re at o dds with the work 
of iDenissenkov fc Pinsonneaui^ l|2008bl ). who found it nec- 
essary to introduce a reduced efficiency of mixing to ex- 
plain the observed abundance trends in CEMP stars. Their 
models produce a ^ inversion in a layer of the star that is 
sufficiently cool that l ittle CN-cycling takes place (see th e 
lower panel in fig. 6 of IDenissenkov fc Pinsonneaultll20()8al ) . 
As such, they find the mechanism is unable to affect the 
surface CN abundances, unless the mixing takes place to 
a greater depth. They postulate that if mixing is able to 
overshoot to around one pressure scale height below the 
inversion, then sufficient CN cycling can take place to af- 
fect the surface abundances. We note that of the previous 
implem ent ations of thermohalin e mixi ng bv lEggleton et ahl 
l|2008ll and ICharbonnel fc ZahnI (|2007l ') did not need to in- 
clude any such overshooting in order to reproduce observed 
abundance variations. 

In Figure|S]we plot the mass fraction profiles of the CNO 
elements and '^He as a function of radius for our carbon- 
normal model. The profile is a very good match to that of 



Denissenkov & Pinsonneault (cf. the lower panel of fig. 6 in 
their paper). We also confirm their finding that the radius at 
whic h the minimum /j is f ound does not vary much with time 
( De nissenkov fc Pinsonne ault 2008bi) ■ However, we find that 
our ^ inversion develops at r/R© = 0.06, which is slightly - 
but significantly - lower than their value of r/R© ~ 0.0675. 
Crucially, this means our mixing is able to proceed to a 
depth where significant CN cycling occurs. We therefore do 
not need to invoke any overshooting beyond the point of 
minimum fj, as Denissenkov & Pinsonneault did. The rea- 
son for this variation in the location of the p inversion is 
unclear. A detailed model comparison would be desirable, 
but insufficient information exists at present to make this 
possible. 

While (5/^-mixing appears to be able to describe the 
trends observed in both carbon-enhanced and carbon- 
normal metal-poor stars, we must add one note of caution. It 
has been pointe d out that rotation could suppress t his mix- 
ing mechanism (|Denissenkov fc Pinsonneaultll2008al ). While 
our code does not include rotational physics, it is not clear 
that simply adding the diffusion coefficients for both pro- 
cesses would give physical meaningful results. 3D hydrody- 
namic calculations of rotating fluid layers that are unstable 
to thermohaline mixing are clearly required. 



4 CONCLUSION 

We have investigated abundance changes in metal-poor gi- 
ants in the context of the J/^-mixing. We find that this 
mechanism accounts for the trends observed, namely that 
carbon-normal stars undergo significant depletion of carbon 
on the upper part of the red giant branch while carbon-rich 
stars do not. The reduced carbon depletion in CEMPs is 
a natural consequence of the proposed mixing mechanism 
and the efficiency of the mixing does not have to be re- 
duced in an ad hoc way. For the carbon-normal model, the 
[N/Fe] abundance remains within the observed spread in 
metal-poor stars and the ^^C/^'^C ratio is reasonably well 
reproduced. The two models in which carbon-rich material 
has been accreted from a companion show substantial in- 
creases in nitrogen but their [N/Fe] values remain within 
the observed range. Their ^^C/^^C ratios are brought down 
to around the observed value by the action of the S^i mixing. 

The above results are presented with the caveat that 
rotation has not been taken into account in these simula- 
tions, and this may be able to supress the mixing. The in- 
teraction between thermohaline mixing and rotation (or any 
additional mixing mechanism) must be taken into account 
in future work and hydrodynamical simulations of such in- 
teractions are clearly warranted. 
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